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Abstract: We show that a hot rotating fluid of relativistic chiral fermions possesses a new 
gapless collective mode associated with coherent propagation of energy density and chiral 
density waves along the axis of rotation. This mode, which we call the Chiral Heat Wave, 
emerges due to a mixed gauge-gravitational anomaly. At hnite density the Chiral Heat 
Wave couples to the Chiral Vortical Wave while in the presence of an external magnetic 
field it mixes with the Chiral Magnetic Wave. The coupling of the Chiral Magnetic and 
Chiral Vortical Waves is also demonstrated. We find that the coupled waves - which are 
coherent fluctuations of the vector, axial and energy currents - have generally different 
velocities compared to the velocities of the individual waves. 
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1 Introduction 

Anomalies in chiral matter lead to unusual transport effects in an impressive variety of 
physical systems and energy scales. The Chiral Magnetic Effect (CME) generates a vector 
(electric) current along an external magnetic field in a chirally imbalanced matter [1-3]. 
The Chiral Separation Effect (CSE) implies the existence of an axial (chiral) current along 
the background magnetic field in dense chiral systems [4, 5] . In a rotating fluid or plasma 
of chiral fermions the Chiral Vortical Effects (CVE) lead to the appearance of both vector 
and axial currents along the axis of rotation [6-10]. Einally, both magnetic held and 
global rotation of the chiral matter produce an energy flux parallel to the corresponding 
axes [11, 12]. All these effects originate either from the chiral anomaly or from the mixed 
gauge-gravitational anomaly. 
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The anomalous transport phenomena are expected to be realized in the quark-gluon 
plasma created in heavy-ion collisions [13-15], in liquid Helium ^He-A [16, 17], in Dirac [18- 
20] and Weyl [21, 22] semimetals, in cold atomic gases [23], in Early Universe [24, 25], in 
neutron stars and supernovae [26-28]. Recent reviews on anomalous transport phenomena 
can be found in Ref. [29-33]. 

In a stationary uniform background the anomalies generate steady time-independent 
currents. However, the chiral matter is also known to support certain types of sound-like 
density waves associated with anomalous transport because the anomalous transport laws 
couple vector and axial charge densities and their currents to each other and force their 
perturbations to be interrelated. As a result, the vector and axial charge densities may 
propagate as a common vector-axial density wave. 

In an external magnetic field the anomalous coupling between the vector and axial 
sectors is given by the CME and CSE. The corresponding gapless excitation is called 
the Chiral Magnetic Wave (CMW) [34, 35]. The CMW propagates along the axis of 
magnetic field with the velocity which depends on the magnetic field strength. This wave, 
which may exist even at zero fermion density, was suggested to reveal itself in heavy-ion 
collisions via the electric quadrupole observables [36, 37]. Certain experimental signatures 
consistent with possible existence of the CMW in quark-gluon plasma were indeed found 
recently [38], although there are arguments suggesting that these signatures may have 
another explanation [39]. 

If a finite-density chiral system is set into rotation then the vector and axial charge 
densities are coupled to each other by the CVE. The associated gapless mode is called the 
Chiral Vortical Wave (CVW) [40]. The CVW propagates along the axis of rotation and 
it may in principle be observed in noncentral heavy ion collisions which create rotating 
quark-gluon plasmas. 

Another gapless excitation, the Chiral Alfven Wave [41], corresponds to sound-like 
oscillations of local velocity of a charged chiral fluid in the presence of an external magnetic 
field. Despite the fluid oscillations are transverse with respect to the magnetic field axis, 
the wave itself propagates along magnetic field lines in a close analogy with the usual Alfven 
modes that exist in various ion plasmas. 

An external electric field may also lead to appearance of a gapless mode due to the so- 
called Chiral Electric Separation Effect (CESE). The CESE generates the axial current in 
a chirally imbalanced medium in the presence of an external electric field. The associated 
sound mode should propagate along the axis of the electric field and may be potentially 
observable in heavy-ion collisions [42]. 

In our paper we demonstrate the existence of a new gapless collective mode associated 
with a coherent propagation of the (thermal) energy wave and the (chiral) axial density 
wave in a globally rotating medium of relativistic chiral fermions. We call this collective 
mode the Chiral Heat Wave (CHW). 

The heat wave is substantially different from the magnetic and vortical waves^. Eor 

^ We discuss three chiral waves with visually similar abbreviations (CMW, CVW, CHW) as well as various 
coupled (CMW and/or CVW and/or CHW) waves. Thus, in order to avoid almost inevitable confusion, 
we often use the terms “magnetic wave”, “vortical wave”, “heat wave”, “magnetic-vortical wave” etc. 
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example, the heat wave may be realized in the absence of a magnetic field background 
unlike the magnetic wave. Moreover, the heat wave may propagate in a zero density 
system unlike the chiral vortical wave. Finally, as we mentioned, the pure Chiral Heat 
Wave couples thermal energy waves with axial density waves, while the Chiral Magnetic 
and and Chiral Vortical Waves - in their original formulation - correspond to a coherent 
propagation of vector and axial density waves in the absence of thermal energy waves. 

The structure of this paper is as follows. In Sect. 2 we review in detail both the Chiral 
Magnetic and the Chiral Vortical Waves. In Sect. 3 we introduce the Chiral Heat Wave 
and discuss its basic properties. As we show in Sect. 4 these waves may mix in different 
combinations and propagate as a common vector-axial-energy density wave. The properties 
(velocity, direction, density content) depend on the concrete physical environment (density, 
temperature, global rotation and magnetic field). A special case of zero-temperature energy 
waves (basically, matter waves) is also discussed. In Sect. 5 we describe a special class of 
non-propagating diffusion modes which we call the Dense Hot Spots (DHS). These long- 
wavelength configurations carry nonzero vector charge and energy density while their axial 
charge is zero. They may only appear in rotating systems subjected to an external magnetic 
field. Finally, Sect. 6 is devoted to a summary of our results and discussion. 


2 Chiral Magnetic Wave and Chiral Vortical Wave 
2.1 Chiral Magnetic Wave 

In this Section we briefly review the Chiral Magnetic Wave (CMW) following Ref. [35]. The 
CMW is a collective gapless excitation in a system (fluid) of massless charged fermions in 
the background of external magnetic field B. The CMW appears due to correlated interplay 
of the CME and the CSE, which describe, respectively, the dissipationless transfer of electric 
charge and chiral charge along the magnetic field [1, 4, 5]: 

jv — o'ycS , (2.1) 

Ja = <^AeB , ( 2 . 2 ) 

The CME (2.1) generates electric (vector) current of the fermions, jy = j along the 
direction of magnetic field, while the CSE (2.2) leads to appearance of the chiral (axial) 
current given, respectively, by a sum and a difference of the right-handed (j'^) 

and left-handed {ji) fermionic currents: 


Jv — Jr + 3l ) 
Ja = Jr ~ Jl ■ 


(2.3) 


The strength of these effects is controlled by the corresponding anomalous transport 
coefficients in Eqs. (2.1) and (2.2): 


ay - 


a A = 


27r2 ’ 

Tv 

27r2 ’ 


(2.4) 
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where ixy is the usual (vector) chemical potential which describes the total density of the 
right-handed and left-handed fermions while ha is the axial (chiral) chemical potential 
which describes the difference in their densities, respectively: 

/iy = 

HA = H5 

The superscripts B in Eq. (2.4) indicate that these transport coefficients correspond to the 
background magnetic field B. 

Let us consider small long-wave perturbations in the vector and axial charge densities, 
pv = jy and PA = 3% respectively. These perturbations, 

5pv{x) = xvv^Hvix) -h xva5ha{x) , (2.6) 

5pa{x) = XAV^Hvix) + xaaSha{x) , (2.7) 

are related to the (local) deviations in the corresponding chemical potentials Spa = Ha — Ha 
via the susceptibilities Xab with a,b = V, A. Here the bar over a quantity indicates a volume 
mean of the corresponding quantity^. 

Equations (2.1)-(2.6) should be supplemented with the conservation laws for the vector 
and axial charges. 


= 2 ^^^ + ’ 
= xil^R - fJ'l) ■ 


(2.5) 


^^,jA = dtPv + djy = 0, (2.8) 

dfiJA = 9tPA + QJa = 0 > (2-9) 

where jy = {pv-,3v) 3 a — {PA-,3a)- Notice that, in general, the axial current is not 

conserved due to the chiral anomaly. However, the conservation law (2.9) is valid in the 

absence of the external electric field which is the case considered in this paper. 

For a totally neutral system of fermions, all chemical potential vanish on average: 
Hy = 0 and p-A = 0- Then we notice that py and pv {PA and ha) are the components of 
true (axial) vectors and therefore the covariance of Eqs. (2.6) and (2.7) under the P parity 
transformation [V -A —V and A ^ A) implies the absence of the off-diagonal terms in 
the susceptibility matrix, XAV = XVA = 0. Since we consider a chirally unbroken phase 
of the system, then the diagonal terms should be the same: xw = XAA = Xi where the 
susceptibility y = xd-^l) is a function of the magnetic field [35]. This statement should 
also be true for py ^ 0 provided that the chiral symmetry is unbroken. Therefore, for the 
chirally symmetric system (with pA = 0) Eqs. (2.6) and (2.7) are simplified: 

Spv{x) = X^Pvix), 
dpAix) = X^Pa{x) . 

From Eqs. (2.1), (2.2) and (2.4) we deduce that a perturbation in electric and chiral 
charges leads to appearance of, respectively, perturbations in chiral and electric currents 

^In order to avoid cluttering of our notations, hereafter we omit the bars over most of the mean quantities. 
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along the magnetic field: 


^3v{x) 


( 2 . 11 ) 

^3a{x) 


( 2 . 12 ) 


while the transverse (with respect to the magnetic field) components of these currents are 
zero. 

Let us now consider the uniform and constant magnetic field directed along the z 
axis, B = Bbz- Without loss of generality one can take eB > 0. Differentiating the 
currents (2.11) and (2.12) over 2 and applying the corresponding conservation laws (2.8) 
and (2.9), 


dt6pv + dz5jv = 0 , dtSpA + = 0 , 


(2.13) 


one gets the following system of linear differential equations which relates the perturbations 
in electric and chiral charge densities: 

eB 

dtdpv(x) + —^dz5pA{x) = 0, (2.14) 

2-n^X 

eB 

dt6pA{x) + A-^dz6pv{x) = 0. (2.15) 

This system can easily be diagonalized by differentiating the first (second) equation with 
respect to z or t (with respect to t or 2;, respectively) and combining them together. One 
obtains two gapless (massless) wave equations, 


{dt - vImw^z) ^Pv{x) = 0 , 
{dt 5 pa{x) = 0 , 


which describe a sound-like propagation of perturbations in electric and chiral charges (and 
their currents) along the axis of magnetic field with the velocity 


eB 

^CMW - ^ . 


(2.17) 


This gapless excitation is the Chiral Magnetic Wave (CMW). 

The CMW is a coupled vector-axial density wave which propagates, qualitatively, as 
follows [35]: 


(i) due to the Chiral Magnetic Effect (2.11) a local fluctuation of a chiral (axial) charge 
density generates a fluctuation in vector (electric) current in the direction of the 
magnetic held; 

(ii) in turn, the huctuation in electric current leads, due to the electric charge conser¬ 
vation (2.13), to a huctuation of the electric charge density next point along the 
magnetic-held axis; 
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(iii) next, the Chiral Separation Effect (2.12) implies that the electric charge density 
generates the axial current which is again directed along the magnetic field axis; 

(iv) finally, due to the conservation of the axial charge (2.13) the fluctuation in the chiral 
current creates a fluctuation in the chiral charge density next point along the direction 
of magnetic field and then the whole processes repeats itself. 

Thus the Chiral Magnetic Wave is a chain-like process which involves the vector (electric) 
and axial (chiral) densities and their currents. 

Summing (subtracting) Eq. (2.15) with (from) Eq. (2.14) and using Eq. (2.3) we get 
two equations which describe propagation of fluctuations with definite chirality: 

{dt T VciAwdz) ^Pl,r = 0 , (2-18) 

where the upper (lower) sign corresponds to the left- (right-) handed fermions. There¬ 
fore the fluctuations in the densities of the left-handed (right-handed) fermions propagate 
opposite to (along) the direction of the magnetic field vector eB. 

The structure of the CMW in terms of the individual density waves is schematically 
illustrated in Eig. 1. The energy component of the CMW will be discussed in Section 4.3. 


6p 



6p 



Figure 1. Structure of the right-handed (the upper plot) and the left-handed (the lower plot) 
Chiral Magnetic Waves in terms of the vector (the green thin solid line), axial (the red thick dashed 
line) and energy (the thick solid orange line) densities for a zero-density (/iy = 0) cold (T = 0) 
nonrotating (H = 0) chiral medium in the magnetic field background B ^ 0. The vector v indicates 
the direction of the wave propagation. 

Before proceeding further we would like to make two important remarks, which concern 
both the Chiral Magnetic Wave discussed in this Section and other chiral waves mentioned 
below. 

First, we notice that similarly to the CME and the CSE, the CMW is emerging at any 
strength of the magnetic field B since in the above derivation no assumption about the 
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value of B was made [35] . However, unlike the CME and the CSE, the CMW is a dissipative 
phenomenon. In the studied long-wave limit the dissipative terms should reveal themselves 
in the dispersion relations in the quadratic order in the longitudinal momentum kz- 

^L,R{kz) = l^VcMwkz - iD\\kl + ... , (2.19) 

where the first term describes the propagation of the CMW with velocity (2.17) while the 
second term is responsible for its dissipation with a longitudinal diffusion constant Dy. In 
Eq. (2.19) the ellipsis stand for higher-order longitudinal terms and for transverse terms. 
Dissipative properties of the CMW were discussed in details in effective hydrodynamics 
and in holographic approaches [43] and in chiral kinetic theory [44]. 

Second, in our paper we consider the Chiral Magnetic Wave and other waves in the 
absence of background flow of the chiral fluid. It is important to stress, however, that the 
fluid velocity may contribute to vector, axial and energy currents and affect the results of 
the paper if they are applied to realistic situations similar to the ones realized in quark-gluon 
plasmas. Realistic implementation of the Chiral Magnetic wave in heavy-ion collisions was 
considered in Ref. [45]. 

2.2 Chiral Vortical Wave 

In this Section we briefly review the Chiral Vortical Wave (CVW) following Ref. [40]. We 
consider a rotating fluid of massless fermions in the absence of magnetic held. A global 
rotation of the huid can be expressed in terms of the vorticity 11= x v, where v is the 
local velocity of the huid how. Similarly to the CME and the CSE, the rotation should 
generate the vector and axial currents of the fermions along the axis of rotation [6-10]: 

= jA = ^A^, ( 2 . 20 ) 


where the associated transport coefficients are a follows: 




k'Vk'A 

TT^ 


^A = 


+ 


My+f^A 

27r^ 


( 2 . 21 ) 


and T is the temperature of the huid. Equations (2.20)-(2.21) describe the Chiral Vor¬ 
tical Effects (CVEs) in the hrst order of the angular velocity fl. Throughout this paper 
we assume that the chiral huid of relativistic fermions rotates slowly so that the linear 
approximation (2.20) is valid. 

According to Eqs. (2.20) and (2.21), at hnite temperature the rotating neutral huid gen¬ 
erates the axial current along the rotation axis = T'^fl/6. However, unlike the transport 
coefficients for the CME and the CSE (2.4), the CVE coefficients depend quadratically on 
the chemical potentials and therefore in the neutral huid an analogue of the CMW cannot 
appear. 

However, if at least one of the chemical potentials is nonzero, then Eqs. (2.20) and 
(2.21) indicate that huctuation(s) 5^Ay on top of the corresponding mean value(s) fiA. ,v / 
0 would couple linearly to the current(s) and may potentially lead to a wavelike excitation 
in a manner of the CMW that was discussed in the previous section. However, /i^ cannot 
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be have a nonzero value in a realistic system in thermal equilibrium. Indeed, due to the 
existence of topological chirality-flipping processes the assumption of a nonzero value of 
the chiral density in thermodynamic equilibrium is not physical. Therefore, in our paper 
we always consider a system of chiral fermions with zero average chiral density {^a = 0). 
The vector density may, however, be nonzero {fj,v 7 ^ 0). 

The CVE in a finite-density chirally-neutral rotating fluid generates a steady axial 
current. 


while the average vector current vanish, jy = 0. Expanding the CVE transport coeffi¬ 
cients (2.21) to the linear order in the fluctuations of the chemical potentials 5fiv and S^a 
at fixed temperature T, one gets the following expression for the current fluctuations (2.20): 


OJy — —^OflA 5 - 


n SpA 


X- - ^^y X o - ^^y^ 
03A — 7— 




(2.23) 


Here we have used Eqs. (2.6) and (2.7) assuming - following the line of arguments for 
case of the CMW - that the parity is unbroken in the fluid. In the above equations the 
susceptibility is, in general, a function x = xipv-, l^^l) of the chemical potential pv and 
the angular frequency H. 

Equations (2.23) demonstrate that in a rotating finite-density fluid, a fluctuation in 
the axial (vector) charge couples to the vector (axial) current exactly in the same manner 
as it happens in Eqs. (2.11) and (2.12) which describe the CMW. Thus, the rotating fluid 
should also support a gapless wave-like excitation similar to the CMW. This excitation 
indeed exists and it is called the Chiral Vortical Wave (CVW) [40]. 

The CVW has the same basic features as the CMW. It propagates along the axis of 
rotation ft with the velocity^. 


/UyH 

VCVW — n ) 


(2.24) 


which can easily be derived from Eq. (2.17) by noticing that Eqs. (2.23) differ from 
Eqs. (2.11) and (2.12) by the simple substitution eB —)• 2pyft. Similarly to the mag¬ 
netic wave, the right-handed (left-handed) chiralities in the vortical wave propagate along 
(opposite to) to the direction of the vorticity vector pyfl. 

The structure of the CVW in terms of the individual vector, axial and energy density 
waves is schematically shown in Eig. 2. The energy density content of the CVW will be 
discussed in Section 4.2. At zero temperature the energy wave in the CVW is, basically, 
the mass wave (4.22) propagating in a uniform finite-density background given by pv 7^ 0- 
The CVW is suggested to exist in the quark-gluon plasma and its signatures are ex¬ 
pected to be found in heavy-ion collisions [40]. 

®Our definition (2.10) of the susceptibility y differs from the one of the susceptibility Xito ol H®!- [^0] by 
the factor of two (y = 2y^o) since the relations Spv/A = X^i^-viA imply &Plir = {xl2)5p,L/R = Xw^IJ-l/r- 
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Figure 2. Structure of the Chiral Vortical Wave in a dense {pv > 0) cold (T = 0) rotating (fi ^ 0) 
chiral medium at vanishing magnetic field {B = 0). If pv < 0 then the energy density and the 
vector charge density have mutually opposite signs while the relative sign of the axial and vector 
densities remains the same. The notations are the same as in Fig. 1. 

3 Chiral Heat Wave 

3.1 Nondissipative energy transfer 

The energy may also be transferred in a nondissipative way due to a mixed gauge-gravitational 
anomaly. The energy flux of a rotating fluid of massless fermions in a magnetic-field back¬ 
ground is [10, 11, 31]: 


Je — crfeS -h cjgfl, 


(3.1) 


where the energy current is given by components of the energy momentum tensor 




(3.2) 


and the anomalous conductivities are expressed via both chemical potentials and temper¬ 
ature: 


^>3 




PVPA ; 




PA 

3 


^ {Spy -h p\) + 


(3.3) 

(3.4) 


The nondissipative energy transfer may only take place in the presence of a chiral imbalance 
(i.e., with PA 7 ^ 0) because the energy current is a vector while both vorticity and 
the magnetic field B are pseudovectors. Thus, the energy current should only be related 
to n and S by a coefficient which is linear in the pseudoscalar chemical potential pA- 
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(3.5) 


The energy conservation implies 

= dte + djj, = 0 , 

where 

e = , (3.6) 

is the (thermal) energy density. 


3.2 Pure Chiral Heat Wave 


The aim of this paper is to find a new gapless wave-like excitation related to the anomalous 
energy transfer. In our context the energy waves should necessarily be related to local 
temperature fluctuations. Let us consider first a small deviation of temperature 5T (with 
6T <C T) from its equilibrium value T. In a linear approximation the energy fluctuation 
6e = e{T + 6T) - e{T) is 


5e = cy(5r, 


(3.7) 


where cy = cv{T) is the specific heat capacity: 


cy = 



(3.8) 


In order to determine the nature of this wave, we notice that - similarly to the cases of 
the CMW and the CVW - the chiral chemical potential in equilibrium is zero, JIa = 0, so 
that a nonzero value of ^a may be due to fluctuations only. Let us also assume for a moment 
that the usual chemical potential is also zero, /2y = 0 similarly to the simplest case of the 
CMW. Then we find that in the linear approximation the coupling of the dissipationless 
energy current to the magnetic field (3.3) is quadratic in fluctuations so that the vector 
current j a does not enter a linear wave equation that we search for. However, the coupling 
of the energy flow to the vorticity (3.4) has a linear term, Ug = T^5 ^a/ 3, which is sensitive 
to the axial charge fluctuations. Thus, in a finite-temperature rotating fluid the energy 
fluctuations should couple to the chiral charge fluctuations. 

The fluctuations of the chiral charge may also couple to electric charge fluctuations in 
rotating fluid at finite magnetic field via the CMW [35] and in rotating fluid at nonzero 
chemical potential via CVW [40] . In order to demonstrate the existence of the new, energy- 
chiral charge wave, we consider first the rotating (H ^ 0) finite-temperature (T ^ 0) system 
in the absence of magnetic field, B = 0 (thus the CMW does not exist) and at zero density, 
/xy = 0 (thus the CVW does not exist either). In this environment only a pure energy-chiral 
charge wave may propagate as a collective excitation. 

According to the full system of equations (2.20), (2.21) and (3.1), the fluctuations obey 
the following relations: 


Sjv = 

0, 

(3.9) 


TQ 

(3.10) 

^3 A = 

IT*"'- 



(3.11) 

= 

3 


10 



We notice first that the new wave does not have the vector density component (3.9). Then, 
using Eqs. (2.10) and (3.7) we rewrite the system (3.10) and (3.11) in the following form: 



(3.12) 


^3 E — ~T. ^PA 

3x 




which clearly demonstrates the coupling of the chiral (energy) current to the energy (chiral) 


density perturbations, similarly to the coupling of the usual and chiral charge currents and 
their densities by the CMW in Eqs. (2.11) and (2.12), and by the CVW in Eq. (2.23). 

Equations (3.12) indicate that, as the wave moves, 

(i) a local perturbation in the thermal energy 5e is converted into a perturbation in the 
chiral current 

(ii) the perturbation in the chiral current 5jj^ leads to an excess in the chiral charge 
density 5pA the next point along the vector 

(hi) the perturbation in the chiral charge density 5pA in turn, induces a perturbation in 
the energy current 

(iv) the energy current 8j ^ generates immediately a fluctuation in the thermal energy he 
next point etc. 

The process is very similar to the mutual cyclic conversion of the charge density and chiral 
charge density waves in the CMW, apart from the fact that in our case the role of the 
charge density wave is played by the thermal energy (or, heat) wave. Thus, we call this 
new gapless excitation the Chiral Heat Wave (CHW). 

The fluctuations of energy and chiral currents (3.12) propagate along the vorticity 
vector fi. Therefore we set for convenience ft = take > 0 for the sake of definiteness, 
and consider the currents directed along the z axis only. The conservation of the axial 
charge (2.13) and energy (3.5) give us the relations 


dtSpA + d^5j\ = 0 , 
dt8e + dz8j% = 0 , 


(3.13) 


which can now be combined with Eq. (3.12) 


dt8pA + X— dz5e = 0 , 
3cy 


(3.14) 


dt8e + ——dzSpA = 0 . 

3x 




The system of equations (3.14) describes a gapless propagation of the coupled energy and 
chiral density perturbations along the axis of rotation. By combining these first-order 
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equations one gets that both the chiral charge density pA and the energy density pE ^ e 
obey the second-order equations: 


{dt - vluwdz) SpAit, 2 ;) = 0 , 
{dt - ^PE{t, z) = 0, 


where 


VcHW 


I r3 n 
cvX 3 


(3.16) 


is the velocity of the Chiral Heat Wave. Notice that in our derivation we have implicitly 
assumed that the usual thermal diffusion and the axial charge relaxation are so slow so 
that the wave propagates adiabatically. 

Let us now discuss the structure of the CHW in terms of energy and chiral charge 
densities. According to Eq. (3.9) the usual charge density, pv = Pr+Pl, does not propagate 
in this wave so that in the CHW the densities of the right-handed and left-handed fermions 
are always opposite to each other: pR = —pL- A diagonalization of Eqs. (3.14) indicates 
that the linear combinations of the energy density and the chiral density, 


SS±{t, z) = 6e{t, z) ± 



obey, respectively, the following linear equations: 


(3.17) 


{dt i Vctitfidz) 6E±{t, z) — 0 . 


(3.18) 


The combinations £1+ and £- represent the pure Chiral Heat Waves, in which the 
energy and chiral charge densities have, respectively, the same and opposite mutual signs: 


6e{t, z) = ± 



6pA{t,z ), 


(3.19) 


and which propagate, correspondingly, along and against the direction of the vorticity 
vector. As we have mentioned, the vector charge density component in the pure heat wave 
is always zero. The structure of the Chiral Heat Wave is illustrated in Eig. 3. 

According to Eq. (3.19), as temperature T becomes higher the fraction of the thermal 
energy component increases compared to the chiral component. Any combination of (ie 
and 6pA perturbations can be expanded into the individual Chiral Heat Waves 6S± which 
are propagating in opposite directions. 

The dispersion relations for the CHWs in the long-wavelength limit reads as follows: 


^CHw(^2:) 4 i^CHW^2: iDh^ “b 


(3.20) 


where upper and lower signs correspond to £± modes (3.19). Here we have also included 
the dissipative term which was not captured by our linearized analysis (higher-order terms 
are denoted by the ellipsis). 
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Figure 3. Structure of the Chiral Heat Wave in a zero-density {pv = 0) hot (T > 0) rotating 
(r2 7 ^ 0) chiral medium at vanishing magnetic field {B = 0). The notations follow Fig. 1. 


In conclusion of this Section we notice that the chiral heat wave bear (due to its thermal 
nature) a very distant similarity to the “temperature wave” (the “second sound”) of the 
superfluid Helium. We would like also to notice that all three discussed gapless modes 
propagate, in general, with different velocities given by Eqs. (2.17), (2.24) and (3.16) for 
magnetic, vortical and heat waves, respectively. 


4 Mixed waves 


4.1 Chiral Magnetic-Vortical wave 

In this section we point out that the CMW and the CVW should naturally mix with each 
other and constitute a common wave provided the conditions necessary for simultaneous 
existence of both these waves (S / 0, / 0 and pv / 0), are satisfied. In this Section 

we consider the cold medium (T = 0) in order to exclude the presence of the Chiral Heat 
Wave. 

In the magnetic field background B the vector and axial currents of the rotating fluid 
are given by the linear combination of the CME (2.1), the CSE (2.2) and the CVE (2.20): 

jy = GyeB + Oy^ ^ (4-1) 

3A = <yleB + aln . (4.2) 


Since the corresponding relations between the fluctuations of the charge densities and 
currents are linear, we combine Eqs. (2.11) and (2.12) with (2.23) and get 


1 eB^^ 

1 eB^^ 


(4.3) 
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which can again be reduced to Eqs. (2.11) and (2.12) with the substitution B — B^^ 
where the effective magnetic field is: 

eB^^ = eB + 2^^^ . (4.4) 


From Eq. (4.3) we see that the CMW and CVW couple to each other and constitute 
one common wave which propagates with the velocity 


_ 

2'k‘^X 


\eB + 2/iyr2| 
27r2x 


V^CMW + ^cvw + 2 ucmw?^cvw cos (p{B, Q), 


(4.5) 


along the axis given by the direction of the vector B^^. The velocity of the common wave 
depends on the angle (p{B, Cl) between the magnetic field and the axis of rotation. The 
subscript “MV” in Eq. (2.24) refers to the mixing of the Chiral Magnetic and Vortical 
Waves. 

Thus, the mix of the magnetic and vortical waves leads to the change of the propagation 
direction of the common mixed wave as compared to the directions of the individual waves. 
Indeed, the individual CMW wave in the absence of the CVW would propagate along the 
magnetic field B while the individual CVW in the absence of the CMW would propagate 
along the rotation velocity Cl. In the a rotating dense system subjected to magnetic field 
the CMW and CVW always mix and form a common wave which propagates only along 
the common vector (4.4). 

The velocity of the common wave is also changed (2.24) compared to the velocities 
of the pure CMW (2.17) and CVW (2.24). The charge-density structure of the waves 
remain, however, the same: the right-handed modes (with the same vector and axial den¬ 
sities) propagate along the effective magnetic field (4.4) while the left-handed modes (with 
opposite vector and axial densities) propagate opposite to this field. 

The energy density component of the mixed magnetic-vortical wave can be easily found 
from Eqs. (3.1), (3.3) and (3.4): 


eB^^ 

^3e — 27 ^ 2 ^ ~ 1 


(4.6) 


where we have also used the first relation in Eq. (4.3). Equation (4.6) demonstrates that the 
energy component of the common magnetic-vortical wave is, basically, the a mass “slave” 
wave which is tightly bound to the vector charge density wave. In particular, 5e{t, x) = 
x). Thus, the qualitative structure of the mixed Chiral Magnetic/Vortical Wave is 
the same as the structure of the Chiral Vortical Wave (shown in Fig. 2) with the reservation 
that the direction of the mixed wave propagation is collinear to the axis of the effective 
field (4.4). 

One of the possible environments where both CMW and CVW may realize is a noncen¬ 
tral heavy-ion collision. The created fireball of quark-gluon plasma is, basically, a rotating 
fluid of light fermions in a background of strong magnetic field (the latter is created by the 
ion constituents and by the products of their collision). In this case the vectors Cl and B 
are co-aligned with each other and the velocity of the common Magnetic-Vortical wave is 
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a sum of the velocities of the individual CMW and CVW (4.5): 


Vmv — 'I’cMw i Vcvw for S||ri, 


(4.7) 


where the upper (lower) sign corresponds to parallel (antiparallel) orientation of eB and 
Notice that for the specific relation between the magnetic field and the angular velocity, 




(4.8) 


the chiral waves do not propagate in the system at all. Indeed, in the considered environ¬ 
ment (S 7 ^ 0, n / 0 and fj-y / 0) the pure CMW and the pure CVW do not exist alone 
so that the vector and axial density waves propagate only in the form of the mixed Chiral 
Magnetic/Vortical Wave. For the strength of magnetic field (4.8) the velocity of the mixed 
wave is zero, so that no wave propagation occurs. 


4.2 Chiral Vortical-Heat Wave 
4.2.1 Three solutions for dispersions 

In the previous section we have shown that the Chiral Heat Wave emerges in a rotating 
finite-temperature fluid at zero chemical potential in the absence of magnetic field. Let us 
lift a bit our restrictions and consider the same fluid but with a small nonzero chemical 
potential, fiy 7^ 0.^ We already know that the Chiral Vortical Wave should appear in this 
environment as a gapless excitation. These two waves should mix with each other because 
the CHW propagates as the coupled energy and axial (chiral) charge density wave, while 
the CVW couples the vector and axial charge densities. Since the CHW and CVW have 
one common axial channel they should inevitably mix at a finite density. 

Following our previous tactics we use the full system of equations (4.1), (4.2) and (3.1) 
along with Eqs. (2.10) and (3.7) to derive relations between the fluctuations of currents J 
and densities = Q: 

5J = nM6P , (4.9) 


where we used the vector notations in the charge, axial charge and energy space: 


J> = (Q,J) 



(4.10) 


with the matrix 


M = 


/ 0 ^ 
My 

tFx 



(4.11) 


^Here we have explicitly included the prefactors e and /ry since they may take negative values. 

®In order to keep our analysis simple, we mostly consider slowly rotating systems at small chemical 
potential (12 <C T, nv ^ T), so that 0(12^) and 0{^y) terms in energy density (3.6) can be neglected. 
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Using the vector form for the conservation laws of vector and axial charges (2.13) and 
that for the energy (3.5), 

= (4.12) 

one gets from Eq. (4.9): 

(ldt + Md,'^6Q{t,z) = 0. (4.13) 

It is convenient to use the plane-wave ansatz in Eq. (4.13), 

Q(t, z) = (4.14) 

It appears that there are three solutions for the dispersion relation: 

‘^vh(^) = ^Vy-akz , (4.15) 

and 


io{k) = 0. (4.16) 

In Eq. (4.15) the subscript “VH” stands for the mixed Chiral Vortical-Heat Wave. 

Below we discuss the linear dispersion laws (4.15) corresponding to the mixing between 
heat and vortical waves. The special (non propagating) zero-frequency solution (4.16) is 
discussed in Section 5.1. 


4.2.2 Propagating mode: vector-axial-energy wave at T 7 ^ 0 

The dispersion (4.15) corresponds to the velocity of the common. Chiral Vortical and Chiral 
Heat, gapless mode: 


VvH = 





1 


/ ^cvw + 



(4.17) 


Thus, in the rotating dense hot fluid the CVW and the CHW form one vortical-heat wave 
propagating with velocity (4.17) along the vorticity vector ft. Notice that the velocity of 
the mixed CVW/CHW excitation (4.17) is higher compared to the individual velocities of 
the pure CVW (2.24) and CHW (3.16). 

The vector, axial and energy densities inside the mixed CVW/CHW are interrelated: 

6pA{t,z) = ±^^5pv'(f, 2 :), (4.18) 

Vcvw 

6e{t, z) = iuv + ^Pv{t, z ), (4.19) 

where the upper and lower signs correspond to those of Eq. (4.15), the common vortical- 
heat wave velocity Uyn is given in Eq. (4.17) and the pure CVW velocity Ucvw velocity can 
be found in Eq. (2.24). Equations (4.18) and (4.19) correspond to the eigenvectors (4.14) 
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Figure 4. Structure of the mixed Chiral Vortical/Heat Wave in a finite-density {pv > 0) hot 
(T > 0) rotating (H ^ 0) chiral medium at vanishing magnetic field {B = 0). The notations and 
the pv < 0 case are described in the caption of Fig. 2. 

of Eq. (4.13) with the dispersion relations (4.15). As expected, the mixed heat/vortical 
wave is composed of all three vector, axial and energy waves (shown in Fig. 4). 

The relation between the vector and axial charges in the mixed wave (4.18) can also 
be rewritten as follows: 

bpA{t, z) = ±J 1 -b ^1 -b bpv{t, z). (4.20) 

Thus, in the mixed wave the axial charge density is always larger than the vector (usual) 
charge density. 

The uj^ modes - which are primarily composed of the right-handed density wave 
with admixture of the left-handed wave - are propagating along with the direction of 
the rotational velocity (vorticity) ft, while the uj~ modes are predominantly left-handed 
densities which propagate opposite to ft. Both these modes are coherent with energy 
density waves which have the same sign as the charge density (4.19) as it is shown in 
Fig. 4. 

4.2.3 Zero temperature: chiral vortical and mass waves 

In the low-temperature limit T —?■ 0 the pure CHW does not propagate (3.16) as its 
velocity vanishes, Vchw{T = 0) = 0. In this limit the u:~^ and uj~ mixed modes (4.15) 
become, respectively, the pure right-handed and pure left-handed modes carried by the 
CVW. Notice that even in this limit - when the Chiral Heat Wave is absent - the energy 
density in the CVW is nonzero due to the presence of the chemical potential pv 7 ^ 0. 
Basically, the Chiral Vortical Wave induces a mass wave in the cold matter (4.19). The 
mass wave is propagating without mass transfer, at least in the linear approximation. 
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Moreover, according to Eq. (4.17), the emerging mass waves are not affecting the velocity of 
propagation of the CVW provided the temperature is sufficiently low, so that the following 
conditions are both satisfied: 


T <C 


) 


T <€. 


TTX ' 


(4.21) 


The charge and energy density content 
Eqs. (4.18) and (4.19): 


inside the CMW can be easily read from 


6pA = ■ 


(4.22) 


We recover the result of Ref. [40] which shows the existence of the pure right-handed 
and left-handed waves which are propagating in opposite directions. In addition, we have 
found that each of these density waves induces the mass wave which is proportional to the 
chemical potential pv- 


VL — ~VCYW ) 

SpL^O, 

SpR = 0 , 

^ 6e = pv^PL , 


vr — +Ucvw ) 

SpL = 0, 

5PR 7 ^ 0 , 

_ 5e = pvSpR. 


(4.23) 


Summarizing, at sufficiently low temperatures (4.21), neither the chirality nor the 
velocity of the CVWs is affected by the presence of the energy wave induced by the mixed 
gauge-gravitational anomaly. However, at higher temperatures the propagating CVW is 
always accompanied by the thermal (“heat”) energy density wave (4.19). The velocity of 
the mixed vortical-heat wave is higher compared to the velocity of the original vortical wave. 
Moreover, the mixed chiral vortical-heat wave is neither right-handed or left-handed (4.18) 
contrary to the pure CVW. 


4.3 Chiral Magnetic-Heat Wave 

Now let us consider the case of rotating (H ^ 0) finite-temperature (T / 0) neutral {pv = 
0) fluid in external magnetic field {B ^ 0). Since the chemical potential is absent then the 
chiral vortical wave does not exist. The propagation of the chiral magnetic and heat waves 
is described by the following relation between the charge and density fluctuations: 


^3 A 

^Je 


eB 

eB 

2-k‘^X 


^PA , 
6pv + 


TO. 

3cv 


6e, 


T^ft 


5pA . 


(4.24) 

(4.25) 

(4.26) 


With the help of the conservation laws (2.8), (2.9) and (3.5), these equations may be cast 
into the system of linear equations: 


(ldt + {M-d)') 6Q{t,z) = 0, 


(4.27) 
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where the charge vector Q is given in Eq. (4.10) and 


M = 


/ « sSf 0 

eB p, rn 

V 0 ^ 0 


The dispersion relation for the density waves in the system reads as follows: 

{Vijk^k^ -uj^)uj = 0, 


(4.28) 


(4.29) 


where the matrix in the momentum space 

ft) = efef + efef, (4.30) 

depends on two unit vectors directed along the magnetic field and the angular velocity, 
respectively: 


e 


B 


B 

~B ’ 



(4.31) 


In Eq. (4.30) the velocities of the pure magnetic and pure heat waves are given in Eqs. (2.17) 
and (3.16), respectively. 

There are three branches of solutions of Eq. (4.29) for the dispersion relations: 


and 



(4.32) 


cj = 0 . (4.33) 

The trivial solution (4.33) is considered in Sect. 5.2 below. 

In order to deal with the nontrivial dispersion relation (4.32) we simplify (without loss 
of generality) our calculations by setting the rotation vector ft along the 2 ;-axis and turning 
the vector of magnetic field B in the xz-plane: 


B = {B sin ip, 0, B cos (p), ft = (0, 0,11), 


(4.34) 


where (p = (p{B, ft) is the angle between the magnetic field B and the rotation axis ft. 
The eigensystem of the matrix (4.30) is was follows: 


Vea = Xaea, a = ±,y, 


(4.35) 


where the trivial eigenvalue Ay = 0 corresponds to the (unit-length) eigenvector By along 
the y axis, which is orthogonal to both ft and B. The nonzero eigenvalues, X± = 
can be expressed via the following quantities: 


2 _ ± I 2 

^MH,± 2 ' ^CMW 






2 

CMW 






. sin^ ip 


(4.36) 
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As we will see below i;MH,± are two principal velocities of the coupled Chiral Magnetic and 
Chiral Heat Waves (denoted by the subscript “MH”). 

The eigenvectors corresponding to the eigenvalues (4.36) are located in the xz plane, 


e+ = (sin 0,0, cos 0), 
e_ = (— cos 6, 0, sin 9 ), 


(4.37) 


where the angle 9 depends on the strength and mutual orientation of the magnetic field 
and the rotation velocity: 


9{B,n) = arctan 



<H,+ (S,n)yl 


tan ip{B, Q) 


(4.38) 


The eigenvectors By and e± form the orthonormal system, BaB), = 5ab with a,b = y. 

It is convenient to re-express the momentum vector k in the orthogonal basis (e_|_, e_. By), 


where 


^ ^ ^ kaf^a ) 

a=±,z 


= kx sin 9 + kz cos 9 , 
k- = —kx cos 9 + kz sin 9 , 


(4.39) 


(4.40) 


and the momentum component ky remains unrotated. Then one gets the following disper¬ 
sion relation from Eq. (4.29) 


u){k) = 


kl + vl,^_k‘i 


where the velocities Umh,± are defined in Eq. (4.36). 
The velocity of the wave propagation is 


v{k) = 


dio{k) 

dk 



OJ 


(4.41) 


(4.42) 


where the matrix V, which operates in the coordinate space, is defined in Eq. (4.30). One 
gets explicitly: 


^L,+fe+e+ + <H,-fe-e- 


(4.43) 


First, we conclude from Eq. (4.43) that for non-collinear Cl and B, the propagation of 
the mixed magnetic/heat wave is not constrained to a singe vector contrary to the cases of 
the pure magnetic wave or pure heat wave. The mixed wave may propagate in the whole 
xz plane spanned on the vectors Cl and B. Second, Eq. (4.43) indicates that the wave 
vector k and the velocity vector v, if even they are constrained to belong to the common 
xz plane, are not parallel unless 


- 20 









(i) The wave vector k is directed along one of the eigenvectors e± (so that either A:+ 
or k- is zero). Then the wave propagates along the vectors e_|_ and e_ with the 
velocities fMH,+ and Umh,- [given in Eq. (4.36)], respectively; 

(ii) One of the fMH,± velocities vanishes (so that sin(^ = 0, implying that the rotation 
axis and the magnetic field are collinear to each other, ri||S); 

(hi ) If r;MH,+ = "yMH,-- This is possible if two conditions are satisfied (4.36); the rotation 
axis and the magnetic field should be perpendicular to each other Cl J- B, so that 
ip = ±vr/2, and the both CMW and CHW velocities should be equal, Ucmw = ^^chw- 

Thus, in general, the direction of the phase velocity, given by the wave vector k and 
the direction of the group velocity, given by the vector v do not coincide with each other, 
which is not unexpected given the anisotropic nature of the medium. In order to illustrate 
this anisotropy, let us consider the the wave vector k directed along the axis of rotation 
ft so that fc||r2||e^. Then from Eqs. (4.36). (4.38) and (4.43) we conclude that the mixed 
wave propagates with the following velocity 

'I’o = vo(sin'd, 0, cos??), (4.44) 


where 


Vo = 


>2 


CHW 


+ v^Mw cos2 v?)2 + sin^ 2ip 


vlnw + cos2 ip 


?? = arctan ^cMw^my^cos^ 

VCHW + ^^CMW COS^ ip 


(4.45) 

(4.46) 


For example, if the angle between the magnetic field and the rotation axis is <y9 = 7r/4 = 
45°, and the parameters of the gas are chosen in such a way that the velocities for the pure 
chiral magnetic (2.17) and heat (3.16) waves are the same, v = Vcmw = Vchw, then the angle 
between the phase and group velocity is ?? ~ 0.32 ~ 18° while the velocity of the mixed 
wave is greater than the velocities of any of its pure constituents: ufin = ^j5/3v 1.3?;. 

We visualize this anisotropic effect in Fig. 5. 


The coupled chiral magnetic-heat wave carries the vector, charge and energy densities 
similarly to the coupled vortex-heat wave. 

Let us consider the case when the magnetic field and angular velocity are collinear 
to each other, S||r2||e2. This case is relevant to noncentral heavy-ion collisions. Setting 
(^ = 0 in Eq. (4.36) we obtain the velocity of the coupled CMW/CHW excitation: 



■\/v 


2 

CMW 


-h V 


2 = 
CHW — 



r302 

9cvx ’ 


(4.47) 


where we have used Eqs. (2.17) and (3.16). In this simple case the mixed wave propagates 
along the common axis of B and Cl, and the directions of the wave and velocity vectors 
are collinear to each other, .B||ri||A:||r;. 
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Figure 5. Illustration of the anisotropic propagation of the mixed Chiral Magnetic/Heat Wave in 
the plane spanned on the vectors of the magnetic held B and the rotation axis Cl. If B and H are 
not collinear, then the direction of the wave vector k and the direction of the wave velocity v may 
not coincide with each other. The angles and d are dehned in Eqs. (4.34) and (4.46), respectively. 


Finally, we mention that in the low-temperature limit the CHW ceases to exist while 
the CMW wave remains unaffected. The ffuctuations in thermal energy decouple from the 
vector and axial ffuctuations according to Eqs. (4.24), (4.25) and (4.26). Thus, the mixed 
magnetic-heat wave becomes the pure Chiral Magnetic Wave at T —)• 0. 


4.4 Chiral Magnetic-Vortical-Heat Wave 


Finally, let us assume that all discussed ingredients are present: we consider hot (T / 0) 
rotating (H / 0) finite-density {^v / 0) fluid subjected to the external magnetic field 
H / 0. The fluctuations of currents 5j^ and the densities 5pa (with a = V,A, E) are now 
related as follows: 




SpA 

2tt'^X 


eB + 


pySpA 


H, 


^Ja 


5pv_ 

2tx‘^X 


eB -|- 


f T5e pvSpv \ 
VScy TT^X ) 


n, 


^3e — 


k-V^PA n , ^PA 

2vrW X 



(4.48) 

(4.49) 

(4.50) 


where we have used Eqs. (2.1), (2.2), (2.3), (2.20), (2.21), (3.1), (3.3) and (3.4). Then the 
conservation laws (2.8), (2.9) and (3.5) give us the familiar relation (4.27) with, however, 
the following new matrix: 


M = 


( ° 

e.B+2/4\/S^ 


V 0 


eSH-2/iy 17 
0 

{eB+2fj.vn)fiv 

2n^X 3x 



(4.51) 


The magnetic field enters Eq. (4.51) only via its linear combination with the angular 
velocity vector (4.4). As we have seen, this is a feature of the CMW/CVW mixing. More¬ 
over, the structure of the matrix (4.51) coincides with the matrix (4.28) which describes 
the mixing of the chiral magnetic and chiral heat waves. 


- 22 - 















Exploring this analogy further, we notice that the excitation spectrum consists of the 
trivial branch cj = 0 and two gapless waves (4.32) with 


V- ■ — 


f 2 £2 I / B £2 £2 , 

/Cj Bj + yt'Mv'J^CHW \&i Cj 


£ 2 ^B£ 2 ^ 


(4.52) 


Here Umv is the velocity of the mixed Chiral Magnetic/Vortex Wave (4.5), is defined in 
Eq. (4.31), is the unit vector in the direction of the effective magnetic held (4.4): 


and 


_ eB + 2fiv^ 
~B^\ ^ \eB + 2^v^\ ’ 


1 ^CHW , 



/ ^ 

2 

[ttTJ 

- 6 \ 

1 Tcv 


(4.53) 


(4.54) 


The case of the trivial dispersion relation (w = 0) is considered in details in Section 5.3. 

Due to the apparent analogy of the full (magnetic/vortex/heat) wave mixing with the 
simpler case of the magnetic/heat wave mixing, the analysis of the full mixing can be easily 
done. Following Section 4.3, we hnd that the eigensystem of the velocity matrix (4.52) 
gives us the principal velocities (squared) Umvh,± = and the principal vectors (4.37) 
determined by the angle 9: 




9{B,n) 


1 / 2 2 

2 ( ^Mv “I" '^cHw T 27 Umv'J2chw cos (3 

('f^Mv + ^CTW + 27 Umv^^chw cos^f - 4(1 - 7^Xv^chw sin^ ^ ^ > 

(1 - 7 ^)^ctw(^) \ VMY{B,il)sunp{B,Q,) 

u^vh^+(S,S7) j 7 Uchw(^ 2 )+ UMv(-B,ri)cos(^(S,n) 


arctan 


1 - 


(4.55) 

(4.56) 


where /3 is an angle between the effective magnetic held B^^, Eq. (4.4), and the angular 
velocity ft, and the factor 7 is given in Eq. (4.54). The mixed waves propagate with 
velocities Umvh,± along the principal vectors (4.37) and (4.56). Here the subscript “MVH” 
stands for the mixing of the Chiral Magnetic-Vortical-Heat Wave. 

All equations of the end of Section 4.3 can now be applied to the full-wave mixing by 
making the substitution B —)> B^^ where the effective held B^^ is given in Eq. (4.4). In 
particular, we conclude that the direction of the wave vector k of the mixed wave and the 
direction of its velocity v do not generally coincide with each other. 

If the vectors B and ft are collinear, then the mixed wave propagates along these 
vectors with (the absolute value of) the velocity 


V 


lli=t 

MVH 


(ucvw i UcMw) T 


1 + 


\ ttT / V Vcvw ^ 


^CHW 


1 

2 


(4.57) 


where upper (lower) sign corresponds to parallel (antiparallel) orientation of the vectors 
eB and //yn, and pure magnetic, vortical and heat velocities are given in Eqs. (2.17), 
(2.24) and (3.16), respectively. 
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Notice that if the magnetic field takes a very specific value (4.8) then the Chiral Mag¬ 
netic and Chiral Vortical Waves disappear completely and the Chiral Heat Wave remains 
the only gapless collective mode in the system. Indeed, at this strength (4.8) the effective 
magnetic field vanishes (4.4), and the matrix M - which determines the propagations 
of fluctuations (4.51) - gets drastically simplified as it has now only two elements corre¬ 
sponding to the Chiral Heat Wave (3.14). Basically, for this value of the magnetic field (4.8), 
the magnetic and vortical waves exactly cancel each other implying Vcvw = =F^^cmw and, 
consequently, Umvh ^ i’chw in Eq. (4.57). The remaining heat wave propagates with the 
standard heat velocity (3.16), which is affected neither by the presence of the magnetic 
field B nor by the rotation ft. 

At the end of this section let us consider briefly the behavior the of mixed magnetic- 
vortical-heat wave in two specific cases. 

The fate of the mixed magnetic-vortical-heat wave in the low-temperature limit can 
easily be seen from the structure of matrix (4.51) which describes the generation of the 
anomalous currents by charge (energy) density fluctuations. At low temperature the heat 
wave disappears while the magnetic and vortical waves still exist and remain coupled to 
each other. The energy current is replaced by the mass current along the direction of the 
effective magnetic field (4.4). The associated mass wave does not, however, influence the 
velocity of the coupled magnetic-vortical wave because in this limit the magnetic-vortical- 
heat velocity (4.55) reduces to the magnetic-vortical velocity (4.5). 

If one keeps temperature nonzero but stops rotation (fi = 0) then both the vortical and 
heat waves formally disappear leaving the CMW alone with the matter wave. According to 
the form of the mixing matrix (4.51), the mass wave is generated by the axial component 
6pA of the CMW. The mass wave influences neither the velocity nor the chiral content of 
the CMW neither at zero temperature nor at finite temperature. 

5 Non-propagating diffusive modes: Dense Hot Spots 

So far we discussed sound-like collective modes which correspond to coherent propagation 
of the vector charge density, axial charge density and/or energy density waves along the 
axis of magnetic field and/or angular velocity vector (or the combinations of the latter 
two). In the lowest, linear order in momentum these waves possess the linear dispersion 
relations, oj = ^vkz, with the corresponding velocities v. However, in certain environments 
we have also found the presence of zero-frequency solutions, a; = 0. In this section we dis¬ 
cuss these solutions and demonstrate that they describe certain non-propagating diffusive 
configurations of energy and vector charge densities. These solutions exist due to interplay 
between axial and mixed gauge-gravitational anomalies in the system. 

5.1 Rotating hot dense fluid in the absence of magnetic held 

In Sect. 4.2 we have demonstrated that the vortical wave mixes with the heat wave in 
a rotating finite-density fluid at finite temperature in the absence of magnetic field. In 
addition, we have observed a new branch of solutions corresponding to the identically zero 
dispersion (4.16) in the linear order of the wave vector k. This new mode is obviously a 
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non-propagating object since its velocity is identically zero in the reference frame defined 
by the chemical potential /iy: 




du}{k) 

dk 


= 0 , 


(5.1) 


We call this diffusive mode “the Dense Hot Spot” (DHS) because this mode has an ex¬ 
cess both in the vector charge density and in the thermal energy density. According to 
Eqs. (4.11) and (4.13) the density fluctuations in DHS are related to each other: 


6e{z) = 


2)Cy^V 

vr^Tx 


5pv{z ), 


5pa{z) = 0 . 


(5.2) 


The fluctuations in energy and the vector charge densities have mutually opposite signs 
while the fluctuation in the axial charge density is identically zero in the DHS. 

The mentioned example of the DHS corresponds to fluctuations in vector and energy 
densities (5.2) which (in the linear order in k) do not generate fluctuations of currents: 


^3v — ^3 A — ^3 E — ^ ■ 


(5.3) 


Therefore, the DHS is not a propagating diffusive mode. We would like to stress that the 
DHS is not a standing wave which could be composed from two counter-propagating pure 
vortex and heat waves. 

Anticipating the inevitable dissipation of the vector charge and thermal diffusion, we 
expect that the dispersion relation of the Dense Hot Spots should contain, to the lowest 
order in momenta, the dissipative term only: 

^Diis{kz) = —-|- ... . (5.4) 

Once this fluctuation is created it would diffuse without propagation. The shape of the 
DHS is not determined by the linear anomalous relations, so that for small densities it can 
be an arbitrary function of the longitudinal coordinate z. 

Thus we arrive to the following convenient definition for the Dense Hot Spot: 

(i) The DHS consists of non-propagating diffusive lumps in the energy and vector charge 
densities which are related to each other; 

(hi) The axial charge density in the DHS is identically zero; 

(ii) The DHS does not generate vector, axial and thermal energy currents via anomalous 
transport laws (5.3). 

5.2 Rotating hot zero-density fluid in magnetic field 

Now let us consider the case of rotating (fi ^ 0) finite-temperature (T / 0) neutral {pv = 
0) fluid subjected to an external magnetic field {B / 0). This setup is relevant to the quark- 
gluon plasma created in noncentral heavy-ion collisions as the plasma is set into rotation 
due to noncentrality of the collision, while the charged nature of the colliding ions exposes 
the plasma to the external magnetic field. Geometrically, both the axis magnetic field B 
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and the the angular velocity vector Cl are co-aligned in noncentral collisions. However, for 
the sake of generality, we will first consider below the case of arbitrary orientation of the 
vectors B and Cl. 

As we have already discussed in Sect. 4.3, in this environment (T / 0, fiy = 0) 
B ^ 0, Cl ^ 0) both chiral magnetic and heat waves emerge and they mix with each 
other. In addition to the coupled chiral magnetic-heat wave we have also found a branch 
of the spectrum with the zero dispersion relation (4.33). This is yet another example of a 
non-propagating diffusive mode, the Dense Hot Spot. 

The charge and energy content of the DHS can be found from Eqs. (3.9), (3.10) and 
(3.11). Since the spot should generate no anomalous currents, we find from Eqs. (3.9) 
and (3.11) that the axial density in the spot is zero, 6pA = 0. Equation (3.10) implies 
the following relation between charge and energy fluctuations valid in the linear order of 
momentum k: 



Figure 6. Illustration of the Dense Hot Spots in rotating hot fluid in a background magnetic field. 
Each plot shows the vector density fluctuation 5pv (z-axis magnitude) and the thermal energy 
fluctuation Se (shown by colors from red/hot to blue/cold) in coordinate space. The axial density is 
identically zero, 6pA = 0. The blue dashed (red solid) arrow points to the direction of the magnetic 
field B (the angular velocity Cl). 

The relation (5.5) between the energy density Cfc and the vector charge density pv,k 
depends on the mutual orientation of the magnetic field B, the angular velocity Cl and the 
wave vector k of the DHS. In Eq. (5.5) the subscript k indicates that the energy density 
fluctuation 8ek,{x) = Cp, cos((A: • x) + , and similarly the vector density 6pv,k, is defined 

for certain wave vector k. Illustrations of the diffusive Dense Hot Spots for a set of mutual 
orientations of the angular velocity and the magnetic field are shown in Fig. 6. 

5.3 Rotating hot dense fluid in magnetic field 

Finally, let us consider the most general situation, when the density, temperature, magnetic 
field and angular velocity are all nonzero. The analysis of the DHS can be done similarly 
to the analysis of the previous section. According to Eqs. (4.48), (4.49) and (4.50) one 
arrives to a density-energy constraint similar to Eq. (5.5), in which the magnetic field B 
is replaced by the effective magnetic field (4.4): eB^^ = eB + 2p,vCl. Thus, we come to a 
conclusion that the vector density fluctuations and the energy density fluctuations in the 
DHSs of rotating hot fluids in the background of magnetic field are qualitatively the same 
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(up to the redefinition of the magnetic field, B —)■ for the cases of zero (/ry = 0) and 

nonzero {pLy / 0) background densities. 

In addition of the DHS-like fluctuations, the rotating dense system of chiral fermions 
may host a “classical” non-propagating DHS mode, in which the magnitudes of vector 
charge and energy densities are not limited to small values. Such “giant” spots are realized 
in a special case, when the magnetic field and angular velocities are collinear to each other: 
S||fi. At zero axial chemical potential, pA = 0, the anomalous vector current (2.1), (2.4), 
(2.20), (2.21) and the anomalous energy current (3.1), (3.3), (3.4) are both zero so that 
the non-propagation condition for these quantities is satisfied automatically. As for the 
anomalous axial current, it is given by Eqs. (2.2), (2.4), (2.20) and (2.21): 

Taking B = Be^ and ft = fle^ along the same axis, one finds that the axial current (5.6) 
vanishes in the DHS if the local temperature and local value of the chemical potential 
satisfy the following relation: 

3pv{x)eB + 3/iy(x)D -|- 7r^T^(x)D = 0 . (5.7) 

Similarly to equations for the chiral waves, this relation for the chiral diffusive spot is 
applicable only in the long-wavelength limit. 

Notice that according to Eq. (5.7) the “giant DHS” at zero temperature is realized at 
the particular value of the chemical potential: 

eB . . 

My = . (5.8) 

Summarizing this section we conclude that the Dense Hot Spots are static diffusive 
modes which appear due to interplay of the axial and gauge-gravitational anomalies in the 
environment that couples the heat wave either to the vortex wave (Section 5.1), or to the 
magnetic wave (Section 5.2) or to the common magnetic-vortex wave (Section 5.3). The 
DHS may only exist in the rotating chiral medium subjected to an external magnetic field. 

6 Summary and conclusions 

We have demonstrated that a hot rotating fluid/plasma of chiral fermions possesses a new 
gapless collective excitation, the Chiral Heat Wave, which is associated with the coherent 
propagation of thermal energy density and chiral charge density waves along the axis of 
rotation. The heat wave propagation is a cyclic process: the mixed gauge-gravitational 
anomaly converts a perturbation in the axial charge density into energy current directed 
along the axis of rotation. Then the energy current heats the chiral medium which generates 
an excess of the axial current along the same axis. Finally, the axial current leads again 
to an excess in the axial charge density further along the rotation axis and the processes 
repeats again. 
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At finite density a rotating fluid supports also the Chiral Vortical Wave [40] which 
mixes with the Chiral Heat Wave at nonzero temperature. Moreover, in the presence 
of an external magnetic field the system may also host the Chiral Magnetic Wave [35] 
which, in a rotating fluid, should couple to the heat wave and which may also mix with 
the vortical wave. Since the mentioned vortical and magnetic waves are propagating due 
to similar cyclic conversion of the vector and axial densities, the mixed vortical-heat and 
magnetic-heat waves involve fluctuations of all three (vector, axial and energy) currents 
and their densities, which appear in different proportions depending on external conditions 
(temperature, density, rotation and magnetic field). 



T / 0 

ft = 0 

ri 7^ 0 

B = 0 

B^O 

B = 0 

B / 0 

PV = 0 

- 

M (2.17) 

H (3.16) 

M+H (4.36), A, DHS 

Pv 

- 

M+p (2.17) 

V+H (4.17), DHS 

M+V+H (4.55), A, DHS 



T = 0 


n = 0 

n 7^ 0 


B = 0 

BAO 

B = 0 

B 7^ 0 

pv = 0 

- 

M (2.17) 

- 

M (2.17) 

pv A 0 

- 

M+p (2.17) 

M+p (2.24) 

M+M+p (4.5), DHS 


Table 1. Physical conditions required for existence of the Chiral Magnetic (M), Vortical (V), Heat 
(H) and Mass (^) Waves, various mixings of these waves and the Dense Hot Spots (DHS) at finite 
temperature T ^ 0 (the upper table) and zero temperature T = 0 (the lower table). The numbers in 
round brackets point out to expressions for the corresponding velocities. The pure Chiral Magnetic 
Wave and Chiral Heat Wave exist at zero density only. Strictly speaking, the Chiral Vortical Wave 
does not exist alone as it is always coupled either to the Chiral Heat Wave at finite temperature 
or to the mass wave at zero temperature. The admixture of the mass (p) wave to any other (M, 
V, H) wave does not influence the velocity of the mixed wave propagation. The label A indicates 
that the Chiral Magnetic-Heat Wave propagates anisotropically as its phase and group velocities 
are pointing out, generally, to different directions. 


In it important to stress that the mixed waves have, in general, different velocities 
compared to the velocities of the individual constituent waves. For example, if the magnetic 
field B and the angular velocity ft are pointing to the same direction (and assuming 
for simplicity that the electric charge e and vector chemical potential pv are both non¬ 
negative), then the velocities of the mixed magnetic-vortical (4.5), magnetic-heat (4.47), 
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vortical-heat (4.17) and magnetic-vortical-heat (4.57) waves are, respectively, as follows: 


where 


'^MV — '*^CMW 4" 'l^CVW ) 


= '/l 


+ vi 


CMW ' ^CHW ) 




+ 






1 + 


I' Hv ^ _j_ ^CMW \ 
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l^cvw J 


^CHW 4~ ("^CMW 4“ '^cvw') 


( 6 . 1 ) 
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2-k‘^X ’ 
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■^^CHW — 


I r3 n 

CVX 3 ’ 


( 6 . 2 ) 


are the velocities of the “pure” Chiral Magnetic (2.17), Chiral Vortical (2.24) and Chiral 
Heat (3.16) Waves, respectively. In Eq. (6.1) the superscript “||” indicates that the ve¬ 
locities are shown for the special case when the magnetic field and the angular velocity 
are parallel to each other, .B||n. We consider slowly rotating system at small chemical 
potential (H <C T, <C T), so that O(H^) and 0(//y) terms in energy density (3.6) are 
neglected. 

Equations (6.1) indicate that the mixing of the Chiral Magnetic Wave with the Chiral 
Heat Wave makes the velocity of the original CMW higher. In other words, the magnetic 
wave (which appears as a result of the axial anomaly) propagates faster at hnite tempera¬ 
ture due to the presence of the mixed gauge-gravitational anomaly. The same is true for the 
mix of the Chiral Vortical Wave and the Heat Wave: the coupling to the energy density 
wave makes the vortical wave faster. Notice that, strictly speaking, the Chiral Vortical 
Wave never exists alone: it is always coupled either to the Chiral Heat Wave (at finite tem¬ 
perature) or to the mass wave (at zero temperature). Due to this inevitable coupling the 
velocity of the Chiral Vortical Wave at finite temperature (4.17) is always higher compared 
to the zero-temperature expression or to idealized formula for a pure Chiral Vortical Wave 
given by the second formula in Eq. (6.2)]. 

As for the mix of the magnetic and vortical waves, the result of the coupling between 
these waves depends on the relative signs of the magnetic field and the angular velocity. 
The mixed wave may propagate faster, slower or even stop propagating at all. The latter 
happens if the the magnetic field and the angular frequency obey Eq. (4.8). There are also 
other effects of the wave mixing. If the angular velocity of the chiral fluid is not collinear to 
the axis of magnetic field, then there exists a mixed Heat-Magnetic wave which propagates 
anisotropically: its phase and group velocities are, in general, not parallel to each other (in 
other words, the wave vector k and the velocity v of the mixed wave are not collinear as it is 
illustrated in Eig. 5 of Section 4.3). The same is true for the triple, Heat-Magnetic-Vortical 
wave mixing discussed in Section 4.4. As for the Heat-Vortical wave (Section 4.2) and the 
Magnetic-Vortical wave (Section 4.1), their phase and group velocities are always parallel 
to each other. 
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A mixing of the heat wave either with the magnetic wave or with the vortical wave or 
with both these waves leads also to appearance of the diffusive modes, the Dense Hot Spots 
which are non-propagating thermal fluctuations with zero chiral charge density but with 
nonzero vector charge density. In the first, linear order in momentum, the DHSs possess 
identically zero dispersion law (w = 0) so that the corresponding phase and group velocities 
are identically zero. The fluctuations in energy density and in vector charge density of a 
DHS are related to each other in such a way that they generate no anomalous vector, axial 
and energy currents. 

In Table 1 we briefly summarize the physical conditions at which the pure Chiral 
Magnetic/Vortical/Heat Waves, the corresponding inter-wave mixings and the diffusive 
Dense Hot Spots may appear. 

We expect that our results should be relevant to noncentral heavy-ion collisions which 
create rotating fireballs of hot quark-gluon plasma subjected to a strong magnetic field. In 
this environment all discussed sound-like modes, namely the Chiral Magnetic, Chiral Vorti¬ 
cal and Chiral Heat Waves should exist. We expect that these waves should inevitably mix 
with each other and form a single collective wave which could be either Chiral Magnetic- 
Heat or Vortical-Heat or Magnetic-Vortical-Heat Wave, depending on the baryon density, 
the strength of magnetic field and the angular velocity of the fireball (c/. Table I). The 
common collective wave propagates as a coherent excitation in vector, axial and thermal 
energy densities. Since in a typical noncentral heavy-ion collision the axis of the magnetic 
field is co-aligned with the angular velocity, the group and phase velocities of the common 
collective wave should coincide. 
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